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ABSTRACT 

Constitutive  equations  of  finite  nonlocal  elasticity  are 
obtained.  Thermodynamic  restriction  are  studied.  The  linear  theory 
is  given  for  anisotropic  and  isotropic  solids.  The  physical  and 
mathematical  properties  of  the  nonlocal  elastic  moduli  are  explored 
through  lattice  dynamics  and  dispersive  wave  propagations.  The 
theory  is  applied  to  the  problems  of  surface  waves,  screw  dislocation 
and  a  crack.  Excellent  agreements  with  the  results  known  in  atomic 
lattice  dynamics  and  experiments  display  the  power  and  potential  of 
the  theory. 


1,  Introduction 


All  physical  theories  possess  certain  domain  of  applicability 
outside  of  which  they  fail  to  predict  the  physical  phenomena  with 
reasonable  accuracy.  While  the  boundaries  of  these  domains  are  not 
known  precisely,  often  the  failure  of  a  given  mathematical  model  is 
indicated  by  its  predictions  that  deviate  considerably  from  experimental 
results  or  dramatically  displayed  by  mathematical  singularities  that  it 
leads  to. 

The  domain  of  applicability  of  a  theory  is  a  function  of  some  internal 
characteristic  length  and  time  scales  of  the  media  for  which  it  is  constructed. 
When  these  scales  are  sufficiently  small  as  compared  to  the  corresponding 
external  scales,  then  the  classical  field  theories  give  successful 
results.  Otherwise,  they  fail. 

Such  is  the  situation  with  the  classical  elasticity  theory 
which  possesses  no  internal  scales.  Yet  all  elastic  materials  possess 
inner  structures  in  the  molecular  and  atomic  scales.  Consequently, 
when  the  external  scales  (such  as  wave  length,  period,  the  size  of  the 
area  over  which  applied  loads  are  continuous),  becomes  comparable  with 
the  inner  scales  (such  as  granular  distance,  relaxation  time,  lattice 
parameter),  the  theory  fails  to  apply. 

In  classical  elasticity,  this  situation  is  demonstrated  dramatically 
by  the  singular  stress  field  predicted  at  a  sharp  crack  tip  and' the  phase 
velocities  that  do  not  depend  on  wave  lengths  of  propagating  waves. 

As  a  result  of  the  former  a  perfectly  sensible  physical  criterion  of 
fracture,  based  on  maximum  stress  hypothesis,  was  replaced  by  various 
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ersatzs  (e.g.  Griffith  energy,  J-integral,  etc).  Clearly,  the  infinite 
stress  is  a  sharp  signal  for  the  failure  of  the  theory  rather  than  the 
failure  of  the  fracture  criterion  which  must  be  based  on  the  physical 
concept  of  cohesive  stress.  Regarding  the  phase  velocity,  at  all  wave 
lengths  from  infinite  to  the  atomic  distances,  we  have  ample  experimental 
measurements  of  dispersion  curves.  Only  at  the  very  large  wave  lengths 
is  there  an  agreement  with  classically  predicted  constant  phase  velocity. 

Waves  having  short  wave  lengths  have  been  observed  to  propagate  with 
much  smnaller  phase  velocities  and  in  fact  they  cease  to  propagate  near 
the  boundaries  of  the  Brillouin  zone. 

The  question  arises:  "Should  we  altogether  abandon  classical 

field  theories  and  appeal  to  atomic  theories  only?"  The  answer  depends 
on  the  characteristic  scale  ratios.  Indeed  if  the  motion  of  each  atom 
in  a  body  is  essential  for  the  description  of  a  physical  phenomenon, 
then  the  lattice  dynamics  is  the  only  answer.  If,  on  the  other  hand, 
the  collective  behavior  of  large  number  of  atoms  is  adequate  for  the 
description,  then  continuum  theory  offers  much  simpler  and  practical 
methodology.  Between  these  two  extremes,  there  lies  a  large  domain  full 
of  rich  physical  phenomena. 

Real  materials  possess  a  very  complicated  inner  structure  full  of 
dislocations  and  impurities.  Moreover  the  force  law  among  the  substructure 
is  not  known.  Consequently,  it  is  virtually  impossible  to  carry  out 
calculations  on  the  basis  of  the  atomic  theories.  Even  if  it  were  possible 
to  accomplish  such  voluminous  and  difficult  computations,  results  would 
be  of  no  practical  value.  All  experimental  probes  possess  some  characteristic 
lengths  so  that  they  can  only  measure  statistical  averages.  Consequently 
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we  need  to  calculate  theoretically  certain  statistical  averages  so  that 
comparison  can  be  made  with  experimental  observations.  Hence  we  are 
back  in  the  domain  of  continuum.  Thus,  continuum  theory  makes  sense 
on  its  own  grounds,  provided  it  is  properly  constructed  to  predict  these 
averages  with  sufficient  accuracy. 

Linear  theory  of  nonlocal  elasticity,  which  has  been  proposed  independently 
by  various  authors  [l]-[6],  incorporates  important  features  of  lattice 
dynamics  and  yet  it  contains  classical  elasticity  in  the  long  wave  length 
limit.  It  is  capable  of  addressing  small  as  well  as  large  scale  phenomena. 

Large  number  of  references  on  the  topic  may  be  found  in  [7]-[9].  Inter¬ 
ested  readers  may  also  consult  [10],  [11]  for  the  nonlocal  fluid  dynamics 
and  [12],  [13]  for  nonlocal  electromagnetic  continua. 

Here  I  present  the  theory  of  nonlocal  nonlinear  elasticity  from 
a  continuum  point  of  view.  (See  also  [14]-[17]).  Constitutive  equations 
are  given  in  Section  3,  where  I  employ  the  global  entropy  inequality 
rather  than  the  local  Clausius-Duhem  inequality  to  place  restrictions 
on  the  constitutive  functionals.  In  Section  4,  I  derive  a  special  class 
of  stress-strain  law  for  the  additive  functionals.  Isotropic  solids 
are  studied  in  Section  5  and  linear  theory  is  presented  in  Section  6. 

In  nonlocal  elasticity,  the  stress  at  a  point  is  regarded  as  a  functional 
of  the  strain  tensor.  For  linear,  homogeneous  solids,  this  introduces 
material  moduli  which  are  functions  of  the  distance.  Physical  and  mathe¬ 
matical  properties  of  these  moduli  are  studied  in  Section  7.  Section  8 
gives  the  field  equations. 

Applications  of  the  linear  theory  begins  with  Section  9  to  wave 
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propagation.  Dispersion  curves  are  obtained  for  the  plane  harmonic 
waves  in  an  infinite  solid  and  for  surface  waves.  Results  are  in  excellent 
agreement  with  the  corresponding  ones  obtained  by  means  of  lattice  dynamics. 
In  Section  10,  I  determine  the  stress  distribution  due  to  a  screw-dis¬ 
location.  Cohesive  stress  that  holds  the  atomic  bonds  together  in  a  perfect 
crystal  is  found  to  coincide  with  the  so-called  theoretical  stress  esti¬ 
mated  on  the  basis  of  atomic  theory  or  experiments.  The  last  section 
(Section  11)  treats  the  crack  tip  problem  for  anti-plane  case  (Mode  III). 
Contrary  to  the  classical  result,  the  crack  tip  stress  vanishes  at  the  tip 
and  possesses  a  finite  maximum  near  the  tip.  The  maximum  stress  hypothesis 
of  fracture  can  now  be  restored.  This  enables  us  to  calculate  the  fracture 
toughness  which  is  shown  to  agree  well  with  experimental  results  on  several 
materials . 

These  few  examples  are  sufficient  to  demonstrate  the  power  and 
potential  of  the  theory.  There  exist  several  other  solutions  in  the  liter¬ 
ature,  dealing  with  dislocations,  cracks,  wave  propagations,  defects,  con¬ 
tinuous  distribution  of  dislocations.  They  also  make  successful  pre¬ 
dictions. 

The  purpose  of  this  lecture  is  to  share  my  enthusiasm  with  you  and 
to  draw  your  attention  to  the  exploration  of  these  new  theories. 
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2,  Balance  Laws 

Just  as  in  classical  field  theories,  the  motion  of  a  material  point 
X  in  a  body  B  with  volume  V  ,  enclosed  by  its  surface  9V  ,  is 
described  by  the  mapping 

(2.1)  x  =  x(X,t)  (2.1) 

where  x  ,  at  time  t  ,  is  the  spatial  image  of  X  ,  in  the  deformed 
configuration  B  having  volume  1/  enclosed  within  its  surface  9U  .  We 
employ  rectangular  coordinates  X^  and  x^  to  denote  the  position  of 
X  and  x  respectively,  and  assume  that 

(2.2)  J  =  det(3xk/9XK)  >  0 
throughout  B  ,  so  that  the  inverse  of  (2.1) 

(2.3)  X  =  X(x,  t)  ,  XeB 

exists  and  is  unique. 

Under  some  mild  assumptions,  local  balance  laws  of  continuum  mechanics 
are  valid  for  the  nonlocal  theory.  Thus,  we  assume  that  the  body  is  made  up 
of  single  nonpolar  species  and  it  is  inert*  Moreover,  nonlocal  gravita¬ 
tional  effects  can  be  neglectedo  Under  these  assumptions,  the  nonlocal 
residuals  may  be  dropped  and  we  have  the  usual  balance  laws 
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(2.4) 

0  *  pvk,k  ■  0  • 

(2.5) 

*ki,k  *  p(ft  - 

(2.6) 

=  rJlk  ’ 

(2.7) 

Vi,k  *  \ 

and  corresponding  jump  conditions  which  we  do  not  list  here  (cf.  Eringen 
[7]  >  [15] » [I7])  •  Here,  p,  v^  ,  t^  ,  f£  ,  e  ,  and  h  are,  respectively, 

the  mass  density,  velocity  vector,  stress  tensor,  body  force  density, 
internal  energy  density,  heat  vector  and  the  energy  source  density.  The 
usual  summation  convention  is  valid  on  repeated  indices  and  a  superposed 
dot  (•)  denotes  the  material  time  rate  and  an  index  following  a  comma 
partial  derivative,  e.g. 

•  _ 

V£  3t  +  V£,mVm  ’  v£,k  ”  3xk 

*k,K  3Xk 

In  contrast  to  classical  (local)  field  theories,  the  entropy 
inequality  is  a  global  statement  covering  the  entire  body.  It  is  expressed 
as 


(2.8) 


r 

[pr,  -  (qk/0)  -  (Ph/6)  ]  dv  >  0 

1/ 


where  n  is  the  entropy  density  and  G  >  0  is  the  absolute  temperature. 
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By  eliminating  h  between  (2.7)  and  (2.8),  we  have 


(2.9) 


f 

|  [-  p(i»ne)  *  tu»til  *  |qkejk]  dv  >  o 

1/ 


where  we  also  introduced  the  Helmholtz'  free  energy 


(2.10) 


i p  =  e  -  0n 


For  some  purpose,  it  is  convenient  to  introduce  the  material  tensors 


TKL  =  J  tk£XK,kXL,&  * 

(2.11) 

CKL  =  Xk,KXk,L  * 


QK  J  qk  XK,k  * 
J  =  PQ/P 


which  enables  us  to  write  Eq.  (2.9)  in  the  form 


(2.12) 


f 

W  t-  P0(*  *  *  k  tkAl  *  F  Qk6,^  dv  i  0 

V 


This  is  the  material  form  of  the  global  entropy  inequality  which  is  valid 
for  the  nonlocal  field  theories.  Clearly,  the  Clausius-Duhem  inequality 
which  constitute  the  integrands  of  (2.9)  and  (2.12)  are  much  too  restricted. 
For  example,  it  is  not  valid  for  the  phase  transformation  where  cer¬ 
tain  thermodynamic  quantities  exhibit  discontinuities  and/or  singularities. 
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3,  Constitutive  Equations 

According  to  the  axiom  of  causality,  [6 , 18]  all  physical  processes  are 
the  result  of  motions  (deformations)  of  all  material  points  of  the  body 
in  the  past  up  to  and  including  the  present  time.  For  thermomechanical 
phenomena,  the  independent  variables  are  therefore 


(3.1)  Y»  =  {x»  ,  0*  ,  0»  } 

~  >  * 

where  a  prime  is  used  to  denote  the  values  of  these  functions  at  (X',t-T'), 
e.g. 

x'  =  x(X',t-x')  ,  9’  =  9(X’,t-T’)  , 

(3.2) 

X'e  B  ,  0  <  t*  <  « 

Constitutive  equations  express  the  functional  dependence  of  the  set 

(3.3)  Z  =  fy,  n  ,  trl,  qk) 
at  (X,t)  on  the  set  (3.1),  e.g., 

(3.4)  =  F[x',9',e»K] 


For  inhomogeneous  materials,  ip  also  depends  on  Xf  explicitly. 
Response  functionals,  such  as  F  ,  must  be  form- in variant  under 
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arbitrary  spatial  translations  and  rotations.  This  implies  that  ip  will 
depend  on  x'  and  x  only  through  | x ’ — x |  .  Since  the  distance  can  be 

expressed  as  a  functional  of  C  ,  it  proves  to  be  convenient  to  replace 
x’  in  (3.4)  by  Cj^  =  CKL(X»,t-x')  ,  i.e. 

(3.5)  4>(X,t)  =  F[C^L,6',e'K,X'] 

This  functional  equation  constitutes  a  fundamental  source  from  which  con¬ 
stitutive  equations  of  large  classes  of  nonlocal  memory-dependent  thermo¬ 
mechanical  materials  (solids  and  fluids)  can  be  derived  (cf.  [6]). 

Here  we  consider  nonlocal  elastic  solids.  In  this  case,  9'  is  uniform 
throughout  B  and  the  memory-dependence  is  not  present.  Consequently, 


(3.6) 


^(X,t)  =  T[Ckl(X'),X’;6] 


where  in  denotes  C^(X’,t)  >  i.e.,  t  is  suppressed  for 

brevity. 

We  assume  that  C  (X')  is  continuously  differentiable  and  it  belongs 
KL  — 

to  a  Hilbert  space  with  an  inner  product  defined  by 


(3.7) 


^l’^H 


H(X'-X)  tr(C’C')  dV* 

v 


where  the  influence  function  H  is  a  positive  decreasing  function  of  its 
argument  such  that  the  integral  in  (3.7)  exists  and 


(3.8) 


H(0) 


1 


In  this  case,  the  space  of  functions  C'  is  a  Hilbert  space  H  with  a 
finite  norm  defined  by 
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(3.9)  ||C’||  =  (C’,C')H 

In  a  Hilbert  space,  any  continuous,  linear,  real-valued  function  f(F) 
has  the  unique  Riez-Frechet  representation 

(3.10)  f (F)  =  (F,C’)h 

There  exists  a  great  variety  of  influence  functions.  As  an  example, 
we  mention 


(3.11) 


H(X)  =  exp(-  a  |X|)  , 


a  >  0 


It  is  now  possible  to  calculate  'p  .  Let 


(3.12) 

then 


P  \p  =  F 
0 


(3.13) 


°o*  2 


sf(CklI4l-D  *  w6 


where  is  the  Frechet  derivative  of  F  with  the  norm  defined  by  (3.9) 

It  is  linear  in  .  Consequently,  it  can  be  expressed  in  the  form 


(3.14) 


6  F  = 


6F 


rC»  dV' 


6  CKL  KL 


V 
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where  the  operator  6(  )/6(  )  represents  Frechet  partial  derivative. 
Substituting  (3.15)  into  (2.12),  we  organize  it  into 


(3.15) 


1  (n  +  ||)6  dV  ♦  i 


[—  T 
l6  KL 


-  2 


V 


wr  <&->*  dv’l  6kl  dv 

KL 


r 


r  r 


7  V,k  dv  * 


r  1  _ >  r  _  i  —  C '  ]  dV '  dV  >  0 

V  KL  6  cCKL  Kl1 


V  “  V  V 

where  a  superposed  asterisk  (*)  indicates  interchange  of  X  and  X'  ,e.g.. 


A(X',X)  =  A(X,X' ) 


It  should  be  observed  that  the  kernel  of  the  last  double  integral  in  (3.15)  is 

skew-symmetric  in  X  and  X’  .  Hence  it  vanishes.  Moreover,  this  in- 

•  • 

equality  is  linear  in  0  ,  and  G  ^  .  For  arbitrary  and  independent 

variations  of  these  quantities  throughout  V  ,  this  inequality  cannot  be 
maintained  unless 


(3.16). 


n  = 


9F 

30  ’ 


(3.17) 


(3.18) 


KL 


=  2 


,6F 


(6C^-)  dV’ 


KL 


qk  -  0 


Note  also  that  we  set  0*  =  0  since  0  is  assumed  to  be  uniform  throughout 
V  ,  The  spatial  forms  of  the  constitutive  equations  follow  from  (2.11) 
which  gives  for  the  stress  tensor 


tk£  Pn  *k,KXil,L  TKL 


(3.19) 
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A.  Additive  Functionals 

For  additive  functions,  in  the  sense  of  Friedman  and  Katz  [19],  a 
representation  theorem  exists,  according  to  which 


(4.1) 


P04> 


G(CKL'  CKl>  *’■  P  dV' 


Since  the  integral  of  F 

X  ,  we  set 

(4.2) 


over  the  volume  V 


★ 

G  =  G 


must  be  symmetric  in  X' 


and 


Consequently, 

(4.3) 


KL 


V 


[  (1G.  ) 
lSc  J 


*  +  dG 


8C 


-]  dV* 


KL 


and  the  Cauchy's  stress,  tensor  is  given  by  (3.19). 

Crystalline  materials  possess  certain  symmetry  regulations.  All 
thirty-two  classes  of  perfect  crystals  can  be  characterized  by  a  subgroup 
of  the  full  group  of  orthogonal  transformations  and  translations 
(Bk)  of  the  material  frame  of  reference,  namely 


(4.4) 


-  SKLXL  + 


det  S  =  +_  1 
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In  the  case  of  amorphous  materials,  and  Bj,  may  be  different  for 

different  species.  Here  we  consider  that  {S^}  and  { BK>  are  constants. 

^  is  invariant  under  the  transformations  (4.4)  and  therefore: 

(4.5)  G  =  G  , 

for  all  members  of  and  {B^}  ,  where 

(4.6)  G  =  G(SC'ST,  SCST  ,  SX’  +  B  ,  SX  +  B) 

For  a  given  symmetry  group,  these  expressions  restricts  the  functional  forms 
of  the  constitutive  equations. 

In  the  case  of  homogeneous  materials,  G  will  not  depend  on  B  ,  so 

that 

(4.7)  G  =  G(C' ,  C,  X'-X) 

For  isotropic  materials,  {s}  is  the  full  group  of  orthogonal  transfor¬ 


mations. 


5,  Isotropic  Solids 

For  isotropic  solids,  {S}  is  the  full  group  of  orthogonal 
transformations  and  Eq.  (4.5)  states  that  G  is  a  function  of  the 
variants  of  C  ,  C*  and  R  =  X’-X  ,  i.e. 


(5.1) 


G(I  ,  I  ;  I.) 

a  a  B 


and  a  function  of  0  ,  where 


ct— 1,2,. ..,6 
3=7,8,.. .,10 


(5.2) 


1 

[4 


tr  C  , 

I2  =trC2, 

II 

tr  C3  , 

tr  C2C’, 

Ic  =  R* CR  , 

I,  = 

r*c2r  , 

tr  CC*  , 

I8  =  tr(C2C»2)  , 

D 

XQ  = 

R  CC*  R 

R*  R 

y 

From  (4.2),  it  follows  that 
(S.3)  G(Ia,  fa;  IB) 


G(I  ,  I  ;  I„)  =  G 

v  a’  a’  3 


since  Ig  =  Ig  . 


Substituting  (5.1)  into  (4.3),  we  obtain 
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(5.4) 


KL 


f 

,  I  r 3G_  .  -  3G_  .  3G 

^  .  l3l  °KL  3I2  KL  3I3  KM  ML 


+  31  (CKMCML  +  CKMCML')  3 1  RKRL 
4  5 


+  31  ^RKCLMRM  +  RMCMKRL^  +  31  CKL 
b  / 


+  (c  r*  r*  +  r  c’  c*  l 
31  LLKM  MN  NL  ^LM  MN  NKJ 

O 


+  %  <Wm  *  Wl  dv' 


The  spatial  expression  of  the  stress  tensor  follows  from  (3.19). 

The  exact  constitutive  equations  (5.4)  may  be  useful  in  dealing 
with  finite  deformation  problems.  In  general,  they  are  too  complicated 
for  practical  applications.  However,  for  certain  simple  geometries  and 
loading  and  in  those  problems  where  the  nonlinear  region  is  localized 
(e.g.,  the  crack  tip),  the  solution  may  be  tractable.  Of  course,  there 
exists  some  important  problems  where  the  nonlinearity  is  the  rule  rather 
than  the  exception,  e.g.  phase  transition,  rubber  elasticity. 
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6,  Linear  Theory 


For  linear  theory,  it  is  useful  to  employ  the  linear  strain  measure 


(6.1) 


EKL  2  ^CKL  '  ~  2  (UK,L  +  UL,K) 


where  U^,  is  the  displacement  vector  in  B  .  A  second  degree  polynomial 
for  G  satisfying  Eq.  (4.2)  has  the  form 


(6.2)  G  =  Z  +  ZklEkl+  2klEkL  +  2  J'KLMNEKLEMN 

+  2  EKLMNEKLEMN  +  2  EKLMNEKLEMN  +  2  ZKLMN’EKLEMN 


1  *2 

where  constitutive  coefficients  E 

KLMN 

and  6  and  they  have  the  symmetry  regulations 


are  functions  of  X1- 


(6.3) 


ZKL  =  ZLK  ’  ZKLMN  =  ZLKMN  =  ZKLNM  ’ 


2  2  2  2 
ZKLMN  =  ZLKMN  =  ZKLNM  =  ZMNKL 


Substituting  (6.2)  into  (4.3),  we  obtain 


(6.4) 


KL 


0  0 

ZKL  +  ZKLMNEMN  + 


V I  p  *  rIV  * 

KLMN  MN 


V 


where 


XX 
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(6.5) 


e„,  = 

1  dV 

KL 

J  KL 

V 

0 

f  72 

KLMN 

J  21  KLMN 

V 

1  ,^1 

KLMN 

2  ZKLMN  +' 

i 

dV' 

n 


from  which  and  (6.3),  we  deduce  the  symmetry  regulations 


jp  =  E*^ 

KL  LK  »  ^KLMN 


E°  =  =  E^ 

LKMN  KLNM  HNKL  ’ 


(6.6) 


ZKLMN 


★ 

ZLKMN  =  ZKLNM  =  ZMNKL 


Constitutive  moduli  E^T  ,  are  functions  of  X  and  6  ,  but 

KL  KLMN 

E'  are  functions  of  X’-X  and  0  . 

KLMN 

From  Eq.  (6.4),  it  is  clear  that  the  nonlocal  effects  are 
represented  by  the  volume  integral,  i.e.  when  ^klmn  =  ®  we  obtain 
the  classical  Hookes  law  of  elasticity. 

One  may  be  tempted  to  incorporate  E^y^  into  the  nonlocal 
moduli  by  writing 


(6.7) 


ZKLMN  ZKLMN  +  ZK 


KLMN 


where  6 (X’-X)  is  the  three-dimensional  Dirac-delta  measure.  With  this, 
Eq.  (6.4)  becomes 


'6-«  tkl  *4* 


f 


rKUIN(r->S-e)  dv«’> 


V 
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In  this  form,  the  continuity  requirement  imposed  on  the  stress  functional 
is  violated.  As  we  shall  see,  this  form  is  useful  for  mathematical  treat¬ 
ment,  but  in  some  instances  (for  compact  operators),  it  leads  to  unbounded 
inverses  which  may  or  may  not  be  acceptable  on  physical  grounds. 

In  Eq.  (6.8),  represents  the  initial  stress  present  in  the 

reference  state,  in  the  absence  of  which  we  set  zj^  =  0  . 

The  material  symmetry  group  {S^}  Place  restrictions  on  Z^^  . 
These  are  of  the  form 


(6,9-)  ZKLMN^’e)  SKPSLQSMRSNT  ZPQRT(~,0:) 

where  R  =  X*  -  X  . 

The  spatial  form  of  the  stress  tensor  is  obtained  by  carrying  (6.8) 
into  (3.19)  and  using 


(6.10) 


EKL  ek£  Xk,KX£,L  *  RKL  =  rkJt  *k,KXA,L 

Xk,K  ~  ^MK  +  EMK  +  ^MK^Mk  ’  p/p0  "  1  "  err 


where  R^  and  r^  are  the  material  and  spatial  rotation  tensors,  re¬ 
spectively.  6^  is  the  Kronecker  delta  when  the  spatial  and  material 
coordinates  are  coincident.  In  terms  material  and  spatial  displacement 


vectors  and  u^ 


we  have 


lKL  -  2  CUK,L  '  UL,KJ  ' 


k£ 


2  (uk,£  “  U£,k')  5 


'k£ 


=  2  (\,£  +  u£,k) 


(6.11) 
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With  these,  Eq.  (3.19)  gives 


(6.12) 


^Jl  ^  err^  uok£  +  °om£.*'ekin  +  rkm^  +  °okm  ^eJlm  +  rJlm^ 


ewi(x')  dv(x,) 


where 


(6.13)  aQ]^  Z'OKL6Kk6L£  ,  °ki,mn  ~  ZKU4N6Kk6LJl6Mm6Nn 

For  homogeneous  and  isotropic  materials,  we  can  show  that,  [7  ],  [13], 


(6.14) 


°ok£  "  °o6kJl  * 


kJlmn 


X»  «v.««  +  V'  (K  &o  +  <$,.A  ) 

k£  mn  km  Jin  kn  Jim 


+  X.  r  “(r  r  6.  „  +  r.  r06  ) 

1  m  n  kJl  k  Jl  iror 

_2 

+  X~  r  (r,  r  6.  +  r.  r  6n  +  r„r  6.  +  r.r  6,  1 

2  ^  k  m  Jin  knJlm  Jlmkn  Jin  knr 


,  -4 

+  X,  r  r,  r„  r  r 
3  k  Jl  ra  n 


where 

(6.15) 


r  =  Vk 


and  Oq  is  a  function  of  x  and  6  but  V,  U*  and  X^  are  functions  of 


f=  j x'  -x |  and  6  ,  i.e. 


(6.16) 


o0(x,6)  ,  X'  =  X'(|x'-x|,6)  ,  y»  =  y»(|x'-x|,6) 


X  =  X  lx'-  x|,  6)  ,  a 

a  a  J  7 


The  appearance  of  XQ  in  (6.14)  indicates  that  the  nonlocal  effects  can 
cause  directional  dependence,  even  in  isotropic  solids.  These  terms  are 
expected  to  be  small  as  compared  to  X'  and  y '  because  of  the  strong 
attenuation  of  the  intermolecular  forces  with  the  distance.  Henceforth, 
we  shall  drop  these  terms  and  alse  assume  that  the  reference  state  is 
stress  free,  so  that 


(6.17) 


k£ 


e  (x* ) 
mn  ~ 


dv(x') 


valid  for  linear  anisotropic  solids,  but  for  isotropic  solids 
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7,  Nonlocal  Elastic  Moduli 


The  nonlocal  elastic  moduli 
certain  physically  acceptable  properties: 


are  required  to  possess 


(i)  From  (6.17),  it  is  clear  that 


(7.1) 


Dimension  (Oj^^)  =  stress/ (length) ' 


This  indicates  that  elastic  properties  of  materials  depend  on  some 
internal  characteristic  lengths  (in  general,  three) ,  a  property  which 
is  lacking  from  the  classical  (local)  theory.  If  e(9)  denotes  this 
internal  characteristic  length  for  isotropic  solids,  then  we  can  put 


(7.2)  A'  =  A (0)  0^  (| x'-xl/e)  ,  y'  =  y(6)  a2(|x'-x|/e) 

where  X  and  y  are  the  classical  Lame  constants ,  and  and 
are  the  attenuation  functions. 

The  internal  characteristic  length  e  can  be  taken  to  be  proportional 
to  the  lattice  parameter,  average  granular  distance,  or  some  other 
internal  length,  depending  on  the  inner  structure  of  the  material  and 
the  accuracy  desired.  We  may  put 

(7.3)  e  =  eQ(0)a 

where  a  is  the  internal  characteristic  length  and  eo(6)  a 
non-dimensional  material  function. 
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In  the  limit  as  e  -*■  0  ,  nonlocal  theory  should  convert  to  classical 
(local)  theory.  Hence 


(7.4) 


lim  (aJ,a2) 
£-*■  0 


=  6 (x'-x) 


It  is  also  clear  that  and  a 2  should  assume  their  maxima  at  x'=x  . 
Similar  limits  are  valid  for  the  anisotropic  moduli  (x' -x>®) 


(ii)  The  nonlocal  moduli  must  be  continuous  and  bounded.  This  is  the 
case, at  least  for  e^O  ,  because  of  the  function  space  to  which 
they  belong.  One  may  also  be  tempted  to  consider  compact  operators 
since  they  transform  bounded  sequences  to  sequences  that  have  con¬ 
vergent  subsequences.  This  however  leads  to  unbounded  inverse  in 
an  infinite  dimensional  space.  In  some  physical  situation,  this 
could  be  acceptable  (e.g.  crack  tip,  where  the  strain  may  be  in¬ 


finite)  .  Contrary  to  o. 


k£ran’ 


in  the  original  Eq.  (6.4) 


^KLMN  ,  and  the  corresponding  spatial  moduli  may  define  a  compact 
operator. 


(iii)  The  group  symmetry  of  (x1  -x,0)  is  identical  to  (6.9)  as  is 

clear  from  (6.13).  For  isotropic  solids,  it  is  expected  that  the 
attenuation  of  intermolecular  attractions  are  the  same  in  all  di¬ 
rections.  Consequently,  we  may  also  take 


(7.5) 

subject  to 

(7.6)  Ci  (0) 


a  =  ou=a(|x'-x|/e)  , 


max  a,  lima  =  6  (x'-x) 

c  -*■  0 


(iv)  All  materials,  at  some  internal  length  scale,  must  be  considered 
discrete.  This  characteristic  of  solids  can  be  introduced  by 
means  of  a  cut-off  length  either  through  the  so-called  quasi- 
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continuum  approach  [1],[4],  or  equivalently  using  a  finite  support 


for  the  Fourier  transform  of  the  nonlocal  moduli  o 


k£mn 


and  a  ,  [4],  [13]. 


The  idea  is  the  same  one  that  is  used  in  lattice  dynamics  and 

information  theory:  Let  f (n)  be  a  scalar  or  tensor  valued  function  of 

a  vector  n  which  takes  discrete  values  at  lattice  points  through 

k  k 

all  displacements  of  the  origin  by  vectors  n=n  e^  (n  and  k 
are  integers  e^:  (e^ ,e^ ,6^) ) .  Let  e  be  the  reciprocal  triad  to 

~k  *  ^ ® • 


(7.7) 


k 

e 


*  P 


Define  a  sampling  function 


(7.8)  6b(x)  =  (2tt) 


-3 


exp(ik*x)  dv(k) 


where  the  domain  of  integration  B  is  a  parallelepiped  kj_<7T}  . 

Suppose  that  the  Fourier  transform  f(k)  of  a  continuous  function  f(x) 

Arf  A# 

can  be  represented  by  a  Fourier  series 


(7.9)  f  (k)  =  vn  l  f  (n)  exp (i  n*k)  , 

n 


k  e  B 


where  Vq  is  the  volume  of  the  cell  B  ,  and  f(k)  =0  outside 


of  B  . 
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From  the  inversion  theorem,  it  follows  that 

(7.10)  f(x)  =  v0  l  f  (n)  fiB(x-n) 

n  ~ 

Since,  according  to  (7.8),  6n(0)  =  v_  and  6D(n)  =  0  for  all  other 

D  U  D 

points,  the  continuous  function  f(x)  takes  the  discrete  values  f(n) 
at  lattice  points.  This  representation  can  be  shown  to  be  unique.  Of 
course,  between  the  lattice  points  f(x)  can  be  quite  arbitrary.  Thus, 
by  means  of  the  sampling  function  6^  ,  we  can  replace  a  function  with 
discrete  values  by  a  continuous  function.  6  plays  the  role  of  the 

D 

Dirac-delta  distribution.  The  following  nonlocal  elastic  moduli 

(,*U)  Ck£mn^'"i,e)  =  Ck£mn(e) 

satisfies  all  the  invariance  requirements,  in  addition  to  being  a  con¬ 
tinuous  function. 

Sampling  functions  for  several  crystals  are  known.  For  example,  for 
a  cubic  crystal,  it  has  the  form 

-3  3  -1 

(7.12)  6g(xf-x)  =  it  II  (x!-x.)  sin[TT(x!-x.  )/a] 

j=l  ^ 

where  a  is  the  lattice  constant •  For  hexagonal,  body  and  face-centered 
cubic  lattices,  see  Kotowski  [201 .  For  the  so-called  Debye  continuum 
where  the  Brillouin  zone  is  considered  to  be  spherical,  £  was  calcu- 
lated  by  Kunin  and  coworkers  [21],  [22]. 

(7.13)  6K(r)  =  (</2Ti2r2)  (^  y  -  cos  <r)  , 

r  =  |  x'-x|  ,  k  =  ir/a 
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where  a  is  the  lattice  parameter.  For  the  two-dimensional  case,  6^ 
has  the  form 


(7.14)  6  (x'-x)  =  (k/2  tt  r)  J  (Kr)  ,  <  =  7T/a 


where  is  the  Bessel's  function. 

We  note  that  SR  and  6  given  above,  have  no  singularity  for 
finite  a  ,  but  they  go  into  Dirac-delta  distribution  when  a->-0  , 
as  they  should. 

The  Fourier  transform  of  6  =  1  in  B  ,  and  vanishes  outside. 

D 

Consequently,  we  have 


(7.15) 


where 

(7.16) 


=  0 


when 

when 


a,  n  =  a,  „  (k)  e  (k) 

kJl  k£ran  ~  mn  ~ 


k.cB 

k.  t  B 
1 


From  this,  by  inversion,  we  have 


(7'17>  5U(!9  * 


where  s]c£mn  is  the  Fourier  transform  of  the  compliance  subject  to 


=  6 


kr^£s 


(7.18) 


ak£mn  Smnrs 
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i.e.,  it  is  the  inverse  matrix  to  •  For  isotropic  materials,  it  is 

given  by 

(7*19)  Sk£mn  =  ^  (6km6£n  +  6kn6£m}  "  2-(3X+2-)  6kA6mn 

For  the  case  (7.5),  it  reads 


(7.20) 


sk£mn 


0 

Sk£mn 


(6)  6(k) 


where  s 


0 


k£mn 
constants  X 


has  the  form  (7.19)  with 
and  v  ,  and 


X 


and  y 


replaced  by  Lame 


(7.21) 


6(k)  =  l/a(k) 


(v)  In  calculations,  the  moduli  c (k)  or  s(k)  in  their  full  generality, 
cause  major  mathematical  complexities.  Often  approximate  expressions 
replacing  them  by  power  series  in  k  are  used.  We  give  here  two  ex¬ 
amples  involving  second  degree  expansions  in  k  for  isotropic  solids 


(7.22) 

(7.23) 


a(k)  =  1  -  e2k2 
B(k)  =  1  +  e2k2 


2 

where  k  = k*k  and  e  given  by  (7.3)  is  independent  of  k  .  We 

note  that  the  linear  term  in  k  is  missing  because  of  the  isotropy  and 
center  of  symmetry. 
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In  physical  space,  these  approximations  are  equivalent  to  constitutive 
equations  of  the  forms 

(7‘24)  °U 

for  (7.22)  and 

(7-25)  "u  - 

for  (7.23),  where  a,  is  the  Hookean  (local)  stress  tensor 


(7.26) 


M 


=  X 


e  6.  „ 
rr  k& 


2y  e 


kZ 


Applying  the  operator 


2  2 
(1-eT) 


to  the  exact  constitutive  equation 


(7.27)  tu 


we  obtain  (7.25)  if 


r 

a(x'-x,6)  ou(x')  dv(x') 

V 


(7.28) 


(1  -  e^V^)a  =  5(x'-x) 


This  indicates  that  in  this  approximation,  a  is  a  Green's  function 
for  this  operator.  This  approach  was  used  extensively  by  Eringen  and  his 
coworkers  (cf.  [23]- [25]). 
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(7.29) 


a 


8(x'-x,e)  tk£(x')  dv(x') 


V 


gives  (7.24)  if 


(7.30) 


(1  +  e2V2) g  =  6  (x' -x) 


Even  though  a  and  6  are  given  approximately  by  (7.22)  and  (7.23). 

We  may  employ  either  (7.28)  or  (7.30)  for  the  determination  of  the 
attenuation  function.  For  example,  when  a  is  determined  from  solving 
(7.28),  we  can  then  use  (7.21)  to  determine  8(x) exactly  instead  of  solving 
(7.30)  which  is  an  approximate  equation  satisfied  by  6  .  In  this  sense, 

(7.27),  with  a  given  by  (7.28),  represents  an  exact  nonlocality.  How¬ 
ever,  a  being  limited  to  a  special  class  of  function  represented  by 
Equation  (7.22),  is  a  limited  or  "short"  nonlocality.  Nevertheless,  with 
this  comppromise,  we  have  achieved  a  major  gain  in  that  the  integro- 
partial  differential  equations  for  the  displacement  field  are  now  reduced 
to  singular  partial  differential  equations  through  the  use  of  (7.24)  or 
(7.25). 

Below,  we  give  solutions  of  (7.28)  for  solids  of  infinite  extents. 
For  these  and  other  types  of  moduli,  see  [5],  [17]  and  [25], 

(i)  One  Dimension: 


(2E)"1  exp(-|x'-xj/e) 


(7.31) 


a(|x'-x|) 
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ii)  Two  Dimensions 

(7.32)  a(|x'-x|)  =  (2tt  O'1  C 1  x ' -x | /e^ 

iii)  Three  Dimensions 

(7.33)  a(l*,-x|)  =  (4tt  £2|x'-x|)  1  exp(-  |x'-x|/e) 

Excluding  (7.31),  these  functions  possess  singularity  at  x'=x  and 

-S#  ^ 

they  go  over  the  Dirac  delta  distribution  when  e  0  .  From  (7.51) 
we  have 


(7.34) 

a 

=  (1  +  £2k2)  exp(ikx’) 

and  since  8  = 1/a  , 

we 

have  formally. 

(7.35) 

8 

=  6(x’-x)  -  £2  6"(x*-x) 

which  upon  substitution  into  (7.29),  gives  (7.25).  Thus  verifying  our 
expectation. 

Singular  kernels  such  as  (7.32)  and  (7.33)  may  lead  singular  a,  . 

K  A/ 

(hence  singular  strain)  fields.  For  niay  be  a  continuous  function 

2 

but  V  t^£  may  be  singular  at  some  point.  The  question  arises  whether 
this  is  allowable  in  a  continuum  theory.  Recall  that  in  the  atomic 
scale  there  exists  repulsive  forces  between  ions  when  they  are  close 
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enough  together.  These  forces  become  infinite  when  x'=x  .  Even  in 
macroscopic  levels  at  a  sharp  crack  tip,  classical  elasticity  leads 
to  infinite  strains.  On  these  physical  grounds,  it  seems  that  we  are 
justified  to  employ  such  singular  kernels.  While  these  kernels  do 
not  possess  an  oscillatory  character  observed  in  crystal  physics,  it 
is  possible  to  achieve  such  oscillations  by  adding  a  fourth  degree  term 
to  the  expansions  (7.22)  or  (7.23),  e.g. 

(7.36)  5=1-  e2k2  +  b  k4 

where  b  is  an  appropriate  constant. 

In  Section  9,  we  shall  see  that  these  kernels  for  the  plane  harmonic 
waves,  lead  to  dispersion  relations,  which  are  nearly  coincident  to 
those  obtained  in  lattice  dynamics,  in  the  entire  Brillouin  zone.  Also, 
cohesive  bond  stresses  in  perfect  crystals  calculated  by  means  of  the 
nonlocal  theory  are  in  good  agreement  with  those  known  in  solid  state 
physics  (Section  10  and  11). 
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8,  Field  Equations 

Upon  substituting  (7.27)  into  (2.5)  we  obtain 


(8.1) 

where 

(8.2) 


a,k  0kA^'^  dv^')  +  P(f£  "  “£)  =  0 


0  0 

°kZ  "k5,mnemn  0k£mn  Um,n 


which  assumes  a  single  attenuation  function  for  anisotropic  solids.  We  put 


3a  ,  , . 

3x^  °W(5  "  ' 


3a 

3x£  ak£ 


daU 


3x* 

k 


H 


in  (8.1)  and  convert  the  first  term  to  a  surface  integral  by  means  of 
the  Green-'Gauss  theorem* 


(8.3) 


r 


a(x'-x)  akj^k  (x')  dv(x') 


“  V  dak  *  0(fi  -  v  ■  0 


1/  31/ 

Using  (8.2)  this  gives  the  field  equation  of  the  nonlocal  elasticity 
for  the  displacement  field 


(8.4) 


f 

“f;'-;5  dV(i,:>  •  j  aV-i>  '1  Vt&'J  ‘‘Vi'’ 


l) 


3U 

+  p(£Jl  "  V  =  0 


Here,  the  surface  integral  represents  the  contributions  of  the  surface  stresses 
(e.g.  surface  tension).  Consequently,  nontocal  thzooy  accounti  {,00.  the  6uA{ace. 
phy&ZcA  06  welt. 

For  the  isotropic  solids  we  replace  ac^£mn  =  °k£mn  (6.17) 
leading  to 


(8.6) 


[(X'  +  2u»)  VV«u'  -  y'V  x  V  x  u']dv'  +  p(f  -  ii)  =  0 


without  the  surface  terms.  If  we  assume  (7.5),  this  gives 


(8.7) 


f 

a(|x'-x|)  [  (A  +  2y)  W*u'  -  yVxVxu']dvT  +  p(f  -  u) 
V 


0 


where  A  and  y  are  the  classical  Lame  constraints. 

For  a  we  may  choose  various  proposals  made  in  section  7. 
If  we  employ  the  approximation  leading  to  (7.25)  we  will  have 


(8*8)  °k£,k  +  (1  -  e2v2)(Pf£  “  P'V  =  0 


or  in  terms  of  the  displacement  field 


(8.9)  (A  +  2y)  VV*  u  -yVxyxu+  (1  -  e2V2)  (pf  -  pu) 


Here  we  see  the  advantage  gained  with  the  particular  choice  of  a 
satisfying  (7.28) . 

The  field  equations  obtained  above  must  be  solved  to  determine 
the  displacement  field  u(x,t)  under  appropriate  boundary  and  initial 
conditions.  Boundary  and  initial  conditions  involving  the  displacement 
and  velocity  fields  are  identical  to  those  of  the  classical  theory. 
Boundary  condition  on  tractions  is  based  on  the  stress  field  t^  >  not 


(8.10) 


tk£  nk  t(n)£ 


where 


)  % 


are  the  prescribed  boundary  fractions. 


9,  Propagation  Waves 


1.  Infinite  Media 


Using  classical  decomposition 


(9.1) 


u  =  V<J>  +  V  x  4* 


Vmip  =  0 


Eq.  (8.9)  with  f  =  0  is  satisfied  if 


(9.2) 

(9.3) 

where 


2  2  2  2  •• 

cx  V  <j>  -  (1  -  e  V  )  4> 

2  2  2  2  •• 

c2  -  (1  -  e  ) ip 


0  , 


(9.4) 


c!  =  t(X  ♦  2y)/p]'5  ,  c2  =  (y/p) 


are  the  classical  phase  velocities  of  irrotational  and  equivoluminal 
waves.  For  plane  harmonic  waves,  Eq.  (9.2)  leads  to  the  dispersion 
relations 


(9.5) 


u/^k  =  (1  +  e"k2)"^ 


For  £=0.39  a  ,  where  a  is  the  lattice  parameter,  the  plot  of  wa/c 
versus  ka  is  shown  in  Fig.  1,  where  the  dispersion  curve  for  the  Bom- 
Karman  model  of  lattice  dynamics  is  also  plotted  for  which 


wa/  c, 


(9.6) 


2  sin(^a/2) 
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The  maximum  deviation  of  this  curve  from  the  nonlocal  result  is  less 
than  6%in  the  entire  Brillouin  zone.  It  should  be  remarked  that  Eq„ 
(9.6)  obeys  the  two  natural  conditions 


(9.7) 


da)/dk 


k  =  0 


dto/  dk 


k  =  TT/a 


0 


but  the  nonlocal  result  (9.5)  obeys  only  the  first  one  of  these.  Accord¬ 
ingly,  the  group  velocity  does  not  vanish  at  the  end  of  the  Brillouin 
zone.  However,  this  situation  can  be  remedied  easily  by  taking 


(9.8) 


a(k) 


2.  2 
e  k 


b  k 


instead  of  (7.22)  where  b  is  a  constant  (cf.  Kunin  [8,11],  p.  38). 

Similarly,  the  dispersion  relations  of  the  equivoluminal  waves 
are  obtained  by  means  of  (9.3). 


(ii)  Surface  Waves 

In  the  plane  x^  =  0  ,  ip  =  ipi^  has  a  single  component  ^(x^^jt) 

so  that 


(9.9) 


u. 


3<t  9^ 

9x^  +  9x2  ’ 


9$ _ 9^_ 

U2  dx2  ~  9x1 


The  general  solution  of  (9.2)  and  (9.3)  relevant  to  surface  waves 
are  of  the  form 


<p  =  A  exp[  -kv.x.,  +  ik(x  -  ct)]  , 

(9.10) 

\p  =  B  exp[  -  kv^x-,  +  ik(x9  -  ct)] 
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provided  vq  is  given  by 


(9.11) 


2  2  222-1 
Va  =  1  -  (c/ca)  [1  -  £  k  (c/ca)  ] 


Using  (9.9)  and  (7.27),  the  stress  field  can  be  calculated.  If  we  set 
Z22  =  t21  =  °  °n  tke  sur^ace  x2  =  0  •  we  tke  dispersion  relations 

[  25  ]: 

(9.12)  [(c1/c2)2  (vj  -  1)  +  2]  Cl  +  V2)  -  4  v^2  =  0 


This  can  be  arranged  into  a  quartic  equation 


(9.13) 


Y(a1Y* 


a7y 


+  a3Y  +  a4) 


0 


where 

Y  =  (c/c2)2  , 

al  =  Ye  +  7  (m+D(£k)2  +  T  (l+4m  "  3m2)  (ek)4  +  m(l-m)  (ek)6  , 
a9  =  -  y  +  y  (2m2  -  m  -  3)(ek)2  -  (1+m  -  2m2) (ek)4  , 

(9.4) 

a3  =  J_m+(2"m  -  m2)(ek)2 

l-2v 

a4  “  1  _  m  ’  m  2 (1-v) 


A  root  of  this  equation  for  £  =0.39  a  and  v  =  0.3  is  plotted  • 
against  ka  in  Fig.  2  where  the  result  of  lattice  dynamic  calculations  of 
Wallis  and  Gazs  [26]  for  KC&  are  also  shown.  If  we  note  that  no  free 
constant  is  available  for  matching,  the  agreement  is  beyond  our  expectations. 

The  nonlocal  theory  also  gives  new  types  of  waves  which  are  not 
present  in  local  theory.  However,  these  waves  disappear  in  the  long  wave 
length  regions  (cf.  Kaliski  and  Rymarz  [27]). 
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It  is  important  to  note  that  the  nonlocal  moduli  a(|x'-x|)  is 
appropriate  only  to  homogeneous  and  isotropic  solids.  Half-space  ceases 
to  be  homogeneous  in  the  vicinity  of  the  surface  x2  =  0,  where  in  a  boundary 
layer  of  a  few  atomic  distances,  the  material  is  inhomogeneous  and  therefore 
a  perturbation  is  necessary  in  a(|x'-x|).  The  present  results  are 
however  valid  so  long  as  the  wave  length  is  much  larger  than  this 
boundary  layer  thickness. 
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10,  Screw  Dislocation 

Consider  a  screw  dislocation  located  at  x^  =  0  of  rectangular 
coordinates  (Fig.  3).  The  displacement  field  has  only  single  .component 

=  w(Xj,X2)  which  satisfies  the  equation 

(10.1)  V2w  =  0 

whose  solution  relevant  to  our  problem  may  be  written  in  plane  polar 
coordinates  (r,0)  as 

(10.2)  w  =  ■£-  6 


where  b  is  the  Burger's  vector.  The  non-zero  components  of  are 

given  by 


9w 


yb 


°3i  =  “Sr  =  -fe  sin6’ 

o  1 


(10.3) 
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3w 

3x9 
o  2 


yb 

2iTr 


cos  0 


We  now  carry  these  into 


(10.4) 


(1  '  £  V  }tk£  °k£ 


and  determine  t.«  by  imposing  the  conditions  that  as  r-*-<»  t,  0=  0  and 


Lk£ 

at  r=0  (crack  tip  surface)  t  =  0  . 


kJ, 

[25  ]  (see  also  [28]) 


This  leads  to 
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(10.5) 


tz6 


Opb 

2ttt 


[1-flC^r/e)]  , 


where  is  the  Bessel's  function.  The  plot  of 

U0.6)  Tz0  =  2tt  e  tzQ/y0b  =  p_1[l-p  Kj  CP) 3  >  P  =  r/£ 


versus  p  is  shown  in  Fig.  4,  where  is  also  plotted  classical  elasticity 
solution.  We  see  that  the  stress  at  the  center  p  =  0  is  singular  for  the 
classical  solution,  but  zero  for  the  nonlocal  solution.  The  maximum  stress 
occurs  at  Pc  =  l.l  and  is  given  by 


(10.7)  t  Q  =  0.3995  yrib/27T  ena 

z0  max  K0  0 


If  we  equate  this  to  the  yield  stress  (theoretical)  t  and  write  h  = 
e0  a/0.3995  ,  this  agrees  with  the  estimate  of  Frenkel,  based  on  an  atomic 
model.  For  e0  =  O.39  ,  b/a  =  1  //2  ,  we  obtain  t^/y  =  0.l2  ,  which  com¬ 
pares  with  the  known  value  0.11  for  A£(fcc)  ,  W,  a-Fe(bcc)  and  0.12 
for  NaC£  ,  MgO  (cf.  [29],  p.  160). 

The  strain  energy  per  unit  length  L  in  x^-direction  is  given  by 

rR  r277 


P0^/L  =  E/L 


_1 

2 


r  dr 


be  be  d8 


substituting  from  (10.5)  and  ez0  =  b/27Tr  ,  we  obtain 


E/L  =  (yb  /8r) [£n(R/2e)  +  Kn(R/e)] 
o  u 


(10.8) 
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where  R  is  the  outer  radius  of  the  solid.  Again  we  see  that  for  finite  R 
the  strain  energy  is  finite.  Eq.  (10.8)  shows  how  the  energy  grows  with 
size. 

The  same  problem  can  be  treated  by  means  of  the  nonlocal  theory 
of  continuous  dislocations  [  30  ] .  For  the  screw  dislocations  for  a  given 
Burger's  vector  b  ,  we  have 


(10.9) 


9e 


32 


3x, 


3e 


31 


3x„ 


b 

2 


The  equilibrium  equations  reduce  to 


(10.10) 


!!».!!«.  o 

3x^  3x2 


For  the  stress  tensor  t^  and  t^  >  we  have  the  constitutive  equations 


(io.ll)  It,,.*  }  = 


y(|x'-x| ){e^1,e^2>  da' 


where  the  anti-plane  strains  are  given  by 


(10.12) 


3w 

e31  =  3x^  * 


3w 

e32  -  3x2 


We  introduce  the  stress  potential  <p  by 


(10.13) 


,  - 

r31  '  "  3x2  ’ 


i-  -  il_ 

32  3Xj 


by  means  of  which  Eq.  (10. 10) is  satisfied.  The  Fourier  transforms  of  (10.9) 
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(10.11)  to  (10.15)  may  be  combined  to  give 


(10.14) 


$  =  y  b/k 


2 


y 

where  k  =  (k^k^) 2  is  the  magnitude  of  the  wave  vector. 

Following  Kunin  [8,11],  if  we  employ  the  Debye  Model,  then  the  inverse 
of  (10.14)  is  found  to  be 


(10.15) 


G(|x*-x| )  b(x')  da(x') 


where  the  integration  is  over  the  area  in 


(10.16) 


G(r)  =  |[ 


f  P 


J:(T)  -  1 


(x^,X2)-plane  and 
dt  -  £n(<R)]  , 


p  =  <r  ,  k  =  TT/a 


For  a  single  screw  at  the  origin  r=0  ,  this  gives 


2 

yb  <  x-  JQ(p)  -  1 

O0.17)  t31  -n  2t  -  - 5 - , 

P 

Alternatively, 

(10.18)  T0  =  27Ttz6/yicb 


2 

yb  k  x  Jn(p)  -  1 

t  _  £ _ L  _ 

z32  2tt  2 

P 

[l-JQ(p)]/p  ,  tzr=0 


tzQ  vanishes  at  the  origin  r  =  0.  It  has  slightly  oscillatory 

behavior  for  r>0  (Fig.  5)  and  it  acquires  its  first  maximum  at  p  =  2.76 
with  an  amplitude 


(!0.19)  tz9  max  ~  °*42^  M0bK/2TT 

This  value  is  nearly  identical  to  (10.7)  where  an  entirely  different 
(singular)  kernel  was  used.  Note  however  that  the  location  of  the  maximum 
differ  considerably. 
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11.  Crack  Tip  Problem 

An  infinite  plate  with  a  line  crack  of  length  2c  subjected  to 
loads  at  infinity  is  known  as  the  Griffith  crack.  It  is  well-known  that 
the  classical  elasticity  solution  of  this  problem  gives  stress  singularity 
at  the  crack  tip.  Because  of  this,  no  fracture  criterion  based  on  maximum 
stress  hypothesis  could  be  employed.  Various  alternative  criteria  (Griffith 
energy  criterion,  J-integral,  etc.)  had  to  be  established  to  circumvent  this 
difficulty.  With  the  development  of  nonlocal  theory,  this  problem  has  been 
resolved  by  Eringen  and  his  coworkers  [  31  -  34],  Here  I  discuss  only 
the  so-called  Mode  III  problem,  i.e.  line  crack  subject  to  anti-plane  shear 
load  (Fig.  6). 

The  classical  elasticity  solution  for  a  line  crack  in  a  plate  under 
anti-plane  loading  at  infinity  is  well  known  [35] 

(11.1)  a  =  a23  -  i  al3  =  oQ  z(z2  -  c2)”'* 

where  is  the  applied  shear,  2c  is  the  crack  length  and  z  =  x^+ix0. 


Employing  (11.1)  in  (7.25),  we  determine  the  nonlocal  stress  field 

t,  -r  /e  i6  /2  -i8./2 

(11.2)  t23  -  i  tl3  =  (7T£/2r1)  2  e  (Cj  e  +  C2  e  )  +  o 

where  (r^,0.)  are  the  polar  coordinates  with  the  origin  at  the  right-hand 
crack  tip(Fig.  6).  The  boundary  condition  at  the  crack  tip  is  calculated 
by  considering  the  tip  as  a  small  circular  cylinder  of  radius  r^  so  that 


42 


in  the  limit  t  vanishes  as  r,  +  0  ,  i.e. 

rz  1  9 

(11.3)  lim  t  =0 

~  rz 
0 

This  gives  C2  =  0  and  C.^  =  -  (c/Tre)  2  oQ  and  we  have  [  34  ] 

■iGl 

(11.4)  tz6  '  1  Sr  '  (t32  -  1  V  e 

0 

=  oQ(c/2r1)'*  {(2r2/cr2)^  exp[i  (-6+  -yO]  -  exp(-  r^e)}  expC-i©^) 

(r^,02)  are  the  plane  polar  coordinates  attached  to  the  left  crack  tip 
and  (r,0)  have  the  origin  x^  =  x2  =  0  , 

Along  the  x^-axis  (0  =  0^  =  02  =  0)  ,  t^  acquires  its  maximum 
near  to  the  crack  tip.  From  (11.4),  we  have 

(U.5)  Te(P)  =  tze/o0  =  («)*  tze/Kin 

=  (2p)‘%  [(1  +  YP)(1  +  *£)"*  ~  e"P] 

where 

I, 

(11.6)  p  =  r1/e  ,  Y  =  e/c  ,  Km  =  (ttc)  2  aQ 


It  is  clear  that  t  .  vanishes  at  the  crack  tip  p  =  0  and  acquires  a 


,-4, 


maximum  at  p  =  p^  =  1.2565  since  Y  <<  1  (Y  £  10  )  . 


t  a  =  an(e/c)”2  [(2  p  )2  +  (2p  )"2] 
z0  max  Cr  1  Kc  cJ  J 


(11.7) 
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T.(p)  given  by  (11.5)  is  plotted  against  p  in  Fig.  7.  From  this  it 

U 

is  clear  that 

(a)  The  stress  field  vanishes  at  the  crack  tip  instead  of  being 
singular  as  predicted  by  classical  elasticity. 

(b)  Fracture  begins  when  t  =  t  ,  where  t  is  the  cohesive 

yield  stress.  z0max  y  7 

From  (11.7),  we  have 

(11.8)  Kc/t  =  3.9278  z2 

y 

where  Kc  =  (ttc)  2  is  the  critical  fracture  toughness.  For  e^  =  0.39, 

K  /t  is  shown  in  Table  1  along  with  its  values  based  on  the  classical 

c  y 

y 

fracture  criterion,  namely  Kc  =  (4 y  2  where  Ys  is  the  surface  energy. 
This  table  also  displays  some  experimental  results.  The  present  results  are 
again  in  good  agreement  with  experimental  observations  of  Ohr  and  Chang  [36], 
even  though  further  considerations  are  necessary  for  the  inhomogeneity  of 
the  material  at  the  core  region. 


M aZwuxt 

CZoa>6i£.<xl 

Fluent 

Expe.fume.nt 

Al  (fee) 

1.11 

0.49 

0.3 

Cu  (fee) 

3.86 

0.47 

0.66 

Fe  (bee) 

1.04 

0.42 

0.23 
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(c)  The  maximum  stress  acts  at  p  =  Pq  not  at  the  tip  of  the 
crack  p  =  0.  Although  pc  is  very  small,  this  implies 
that  fracture  initiates  ahead  of  the  crack  joining  the  tip. 

If  the  inhomogeneity  and  the  presence  of  dislocations  near 
the  crack  tip  is  taken  into  account,  we  expect  pc  to  be¬ 
come  larger  than  atomic  dimensions. 

Interactions  of  dislocations  with  crack  was  treated  by  Eringen 
[34],  For  the  solution  of  the  crack  tip  problems  for  modes  I  and  II,  see 
[32],  [33]  and  for  crack  curving,  see  [24],  Point  defects  and  elastic  in¬ 
teractions  were  explored  in  the  works  of  Gairola  [37]  and  Kunin  [8], 
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FIGURE  1:  Dispersion  Curves  for  the  Present  Nonlocal  Model  and  the  Bom-KSrm3n 

Model  of  Lattice  Dynamics 
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Figure  2:  Dispersion  Relations  for  Raleigh  Surface  Waves 


Figure  3:  SCREW  DISLOCATION 
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FIGURE  5.  SHEAR  STRESS  IN  SCREW  DISLOCATION  (Eq.  10.18) 
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Figure  7:  Nondimensional  Shear  in  Mode  III 
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